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Abstract. Here, we analyze the primordial magnetic field transition between a radiative and a matter-dominated 
universe. The gravitational structure formation affects its evolution and energy spectrum. The structure excitation 
can trigger magnetic field amplification and the steepening of its energy density spectrum. 
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CX) , neous magnetic field. In particular, does the structure appear- 
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1. Introduction spectrum? In this paper we refer to the model of the struc- 



The increasing empirical evidence of magnetic fields at 
different cosmological scales raises a series of questions 
on their origin and their dynamical impact on gravi- 
tational structures. Different models of magnetogenesis 
have been proposed (Berezhiani & Dolgov 2004, Grasso & 
Rubinstein 2000, Hogan 1983) and the possible effects on 
structure evolution have also been discussed (e.g. Kim et 
al. 1996, Wasserman 1978). The role of the feedback pro- 
cess (i.e. structure formation on magnetic field evolution) 
is equally important. It is vital to understand whether the 
structure origin and its subsequent evolution may change 



ture origin, presented in Siemieniec & Woszczyna 2004a 
(hereafter AI). It shows the existence of 'explosively' grow- 
ing structures at decoupling, originated in the stochastic 
field of acoustic-type perturbations in the radiation era. 
The present paper deals primarily with the problem of 
the primordial magnetic field transfer at decoupling and 
its later linear evolution. The considered initial magnetic 
field is assumed to be generated in an earlier magneto- 
genesis process. Its final strength, the amplification rate 
at decoupling and the subsequent energy transfer to large 
and/or small scales depend on the velocity field. The pres- 
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ence of velocity flows at the phase of decoupling is a crucial as a bottom- up phenomenon resulting in an 'inverse cas- 

aspect in the evolution of the primordial magnetic fields, cade' type of magnetic energy transfer? Even in the linear 

The magnetic energy may increase at the expense of the regime it seems possible that the velocity field will am- 

kinctic energy of the flows, depending on their specific na- plify magnetic fields soon after recombination and boost 

ture (e.g. shear). the magnetic energy toward larger scales. 

The structure of the paper is following. In Sect. 2 we 

discuss the assumptions and introduce a set of magneto- 2 1 Basic equations and assumptions 

hydrodynamics equations (MHD). In Sect. 3 the solution Tq gtudy ^ we need tQ congider thfi magnctic evomtion 

of the induction equation is given. We discuss the mag- dcscribed by thc induction cquationj which is a part f 

netic energy density spectrum after decoupling in thc last ^ MRD ^ of equations In some cagcs the symmetries 

allow its separation from the MHD system. An example of 



section. 



P = -Pd (1) 
1 



the evolution of selfgravitating magnetized structure with 

2. Preliminary Statements and assumptions planar symmetry (pancake) recently has been analyzed in 

rr,, . ... . . ,. . A . ,. j... Siemicniec & Woszczyna (2004b MCW). Below, we apply 

I he mam purpose of this paper is the investigation of the v ' 

,.,„,,.. . .. . .. ^ thc full set of pressureless MHD equations (corresponding 

primordial held development at decoupling and directly v 

thereafter. Assuming that magnetic fields were created *° 6ClS ' ^ * n ^^^) 

in the early universe by some particle physics processes 

and that the fields were not correlated with the density 

fluctuations, thc following questions are raised. Can the 

correlation between thc magnetic field and matter den- ^ ~ 3 ^ 47tk/> 2o +2lu + 

sity growth be established at recombination, i.e. if ini- ^\j>^ % ^^ ~ 2^^^ ^ 

o 

tially we have a random distribution of magnetic energy, jji — _|_ w ' jjj dH i , (3) 

j -'3 

then, does it become more organized later. The second one 

is how thc final energy distribution is determined by thc where P> ^> H are the density, expansion rate and the mag- 
initial magnetic spectrum. Thc motivation for these ques- net ic strength of the matter; thc hydrodynamic scalars cr 
tions is that thc velocity which exists at recombination and ^ measure the shear and thc rotation, 
affects the magnetic field through the induction equation. In order to study thc cosmological development of the 
Thus if structure formation develops hierarchically, does primordial magnetic field in the linear regime, we solve 
it mean that the cosmological magnctic field also amplifies the above set of equations, assuming the following: 
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1. The evolution of the magnetic field is investigated 5. The perturbations are rotationless but the shear in- 
in the post-recombination phase of a spatially flat FRW troduced by collapse is important. It accounts for the re- 
universe, suits. For compatibility with our earlier paper (MCW), 

we postulate the symmetry in shear tensor of the form 



2. We analyze the magnetic field perturbations (small 
with respect to density perturbations) in the framework of 
the model (AI) describing the density perturbation transi- 
tion between the radiation and the matter-dominated era. 
Thus, we use here the already-derived quantities describ- 
ing the hydrodynamic evolution and structure formation 
which occurs at the transition. 



<j)W = a&Mj/P. (Referring to MCW, where the fluid 
contraction is one-dimensional, it means that a = | for 



the field component which is not contracted (i.e. parallel 
to the velocity) and a = — | for both components lying in 
the pancake plane. As a result of this symmetry the shear 
scalar can be uniquely expressed by the velocity pertur- 
bation in the induction equation. Its first term is equal to 
3. The transition is instantaneous. In this oversimpli- cr-iJ- 7 = ^-(35? 8j — 8ij)W . Thus, if the perturbation 8$ 
ficd approach the recombination coincides with the decou- is known, then the induction equation can be separated 
pling which is not a realistic assumption (at the same time from the whole MHD system. 

one neglects the viscosity phenomena that are present at g The spectral distributions of the hydrodynamic 
this phase). We believe however that the essential evolu- quantities (velocity and density fields) are determined by 
tionary features do not depend very much on this, i.e. the thfi unknown spcctrum f the primordial density pcrtur- 
proposcd instability mechanism persists. bations P(jfe). The primordial magnetic field spectrum is 

also unknown. Both come from the early universe pro- 



cesses and are supposed to be generated independently. 
In our calculations we analyze the scale-free spectra as 
test forms. 

2 2. Propagation of the linear perturbations 



4. Density and velocity perturbations evolve in both 
epochs. Matching of their values at 'decoupling' is per- 
formed according to the Darmois-Israel joining conditions 
(Darmois 1927, Hawking & Ellis 1973, AI). The pertur- 
bations of the primordial magnetic field SH are however 
induced by the density and velocity field perturbations Sp, 
S'd, and since it is attributed to the structure formation The density and velocity contrasts S — 6 = 4^- cal- 

' J J P(0) ' »(0) 

process, they have their origin at decoupling. The magne- culated in (AI) are 'functionals' of the primordial density 

togenesis mechanism producing the primordial seed field spectrum P{k). In the linear regime they are expressed 

is expected to occur prior to recombination i.e. created by Fourier integrals (cf. eq. (38) i (45) in AI), where the 

before the 'initialization' of structure formation. Fourier coefficients {.4k} are functions of the primordial 
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initial spectrum, |.4k| 2 oc P(k), and the density and ve- have a source oc StfHfo, which determines the evolution- 
locity modes in the matter-dominated era are functions of ary behaviour of SH Z . This term itself depends on both 
conformal time 77 and the wave number k given by eqs. the initial injection seed field Hr \ spectrum and the spec- 
(32), (47) in AI. The linearized hydrodynamic evolution, trum of 8$. While the former may be postulated as a free 
eqs. (1) and (2), after introducing the perturbations, parameter, the spectrum S"& is not only a product of the 



P = P(0) + $P (4) 
•d = (o) + 80 (5) 



H l = H\ 0) + SH 



primordial density fluctuation spectrum but also a result 
of changes occurring at the radiation to dust transition. 

The primary result of the transition model (AI) ap- 
plied here is structure excitation with the features of 
is governed by acoustic shock formation at decoupling. The changes in 

g'p _ —fi^Sp — p^Sfl (6) density and velocity field (due to this transition, see fig. 

2 . 2 and 3 in AI) arc felt by the magnetic field. Thus one 

5& = -4irK,8p- -tf(o)fo?, (7) 

can expect that the induced magnetic field SH should be 
substantially correlated with the structure. 

3. Induction equation 

Let us consider the propagation of the magnetic perturba- 
&(0) = -^(0) - 47tk/? (0) , (9) tions (10). For convenience we use here the conformal time 



where the solutions for contrasts (-^-, 4^-) are given by 

V P(0) ' W(0) ; ° ■' 

eqs. (38), (45) in AI. The equations for background dy- 
namics are 



P(0) = P(0)tf(0)i ( 8 ) 

1 

3 l 



• 1 • 1 ,• f \ c sui V = f a (t) 1 dt, where the scale factor 1 a(ri) has the form 

while the linearized induction equation (3) for on , ac- ' J v ; ' y " 

cording to assumption 5, becomes ^ 21 ) in AI ) a ^ = ^Sg^' M is the constant of 

/ o\ 2 motion in the Friedmann equation and 77s represents the 

(SH*) = (a - - J 8VH\ 0) - -ti^SW (10) 

transition time from a radiation to a matter-dominated 
era, i9 = 3 da ^ dt . Equation (10) in the variable 77 can be 



and the seed field satisfies 



#(0) =-^(0)^ ) =>H {0) «a- 2 (11) replaced by 

It is worth noticing from the above that in the case of no — ^ S \ sH 1 )' = — (12) 



asymmetric contractions i.e. /? = a — | = 0, the induc- 



tion equation leads to the known result (identical to the 
background field H/ \ evolution), 5H l oc a~ 2 . Contrary to 1 This form of a(t) is a result of the joining conditions at the 



To obtain its solution in the linear regime we express the 
primordial and perturbed magnetic field as well as the 



that, if P 7^ 0, (for contraction (3 < 0) then in eq. (10) we time of transition with the radiation scale factor 0(77) = y/ 4p ??. 
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velocity perturbation Si} as Fourier series: 



H = r /s 2 (l + r)- 4 / dke ikx /(fc) + c.c. 
SH = n s 2 J dke tkx G(fc,T) +c.c. 

and according to eq. (45) in (AI) 



= ^ = / dke lfex 0(r, k) + c.c. 



Equation (|12[) becomes now 



ur r 



where r = (1 + — ) and 



(13) 
(14) 



(15) 



evolution is determined by the primordial density fluctua- 
tion spectrum. The expansion contrast is defined in the 
dust era by the joining conditions at rys and by the input 
parameter P(k), i.e. the energy spectrum of the acous- 
tic field. According to AI, the general solution for 9 is 
expressed by eq. (45) 



(19) 



(16) 



24tt/3 



where the Fourier coefficients ylk^k x an( i the 

modes Uk(i?) after transition are given by (eq. 47 in AI) 
31/4 

4 6ifc£\ 



+ OC 



2k 
ikYi 
"io"' 



ikYi 



(20) 



dfc'fc' 2 (/(fc- k')e(k',T) + f(k + k')Q*(k',T)). (17) 



where £ = ?y s /V3 . The comparison of Fourier expansion 
The nonhomogeneous term g(k,r) in the induction equa- ( 15 ) with the above allows one to obtain the trans f rm 
tions (16) shows the interaction of the mode H m of the k y For technical rea sons we take in 0(r, k) only the 

primordial field with the velocity modes Q k . This term leading term K T 2 (we drop r -3^ which is equivalent to 
vanishes when = 0. In general, the nonvanishing source the aS y mptotic form of thc so i uti on i.e. to a late time after 



of eq. (16) gives the solution in the form 



recombination. For 0(r, k) one finally gets 



G(k,r)=(^y i C(k) + (^y i £ g(k,r)Q 4 dT, (18) Q(fc, r) = i C x kK 2 e 
where C(k) is to be found from the boundary condition at where C\ = C 3 J S 

2^40 



-ifeS 



(21) 



decoupling i.e. 8H(t = 2) = 0. It is seen from the above 
that the source modifies both the temporal and spatial 
(spectral) evolution of the primordial field. 



The constant C above comes from the primordial fluc- 
tuations. Both quantities forming the convolution in mag- 
netic perturbation source function l|17(l are two 'input pa- 
rameters' depending on the primordial density fluctuation 
The velocity field represented in this calculation by is spectrum P(k) and the spectrum of the primordial mag- 



3 1. Comments on spectral forms 



present at the decoupling time rj = 77s. Together with the netic fields. Both are unknown. The postulates concern- 
density contrast it propagates from the radiation era. Its ing the first quantity commonly refer to the inflationary 
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paradigm. Here we employ no concept for the primordial where 

fluctuation generation. For the sake of mathematical sim- „ , „ „ 

T = 2AttACi(3, 

plicity we adopt the primordial density spectrum of the 

g(k) = il n (k), 

white-noise type i.e. P(k) = const which means Ak = C. fc max 

I n (k) = { dk'((k- fc') 1 - 0' ~ k)^)k' n exp{i(k' - fc)E). 
For the initial primordial field, we do not indicate any ^ 

particular magnetogenesis process with its specific initial Finally, taking into consideration the initial condition 
spectrum . Below we assume that the spectrum generated SH(t = 2) one obtains the solution for G(fc, r) 
by potentially different models of the early magnetogene- ^ _ }_y g(k)(r^ 2 — 4r~ 4 ) (24) 

sis can be described by the power law 

It is important to note the structure of this solution. The 
tih\ — Aun h ^ h ^ h (oo\ temporal dependence of magnetic perturbations includes 

J(K) — A-K K m in < K < K ma x \^^) 

the typical solution of the homogeneous equation (G oc 

f(k) = elsewhere 

r -4 oc a~ 2 ) and the additional component oc r~ 2 oc a . 
In the following calculations we restrict ourselves to the Thus the sheared components, (3 ^ 0, of the induced mag- 
analysis of two cases i.e. n = and n = 1. In the first netic fields are amplified in the linear regime much faster 
case we use the flat spectrum. It is plausible that mag- than those which are not contracted. The spatial depen- 
netic fields have similar strength at different scales when dence is determined by I n (k) i.e. by the convolution of a 
the initial conditions for them were the same at the time primordial magnetic spectrum with the kinematic quan- 
of magnetogenesis. Thus, similar to density fluctuations, tity 9(r, k). 

the white-noise spectrum is suitable. This one and the 

.... ,-, . 4. Magnetic energy density spectrum 

second putative spectrum, j(k) oc k , allow us to put the 

lower spectral limit at k m i n = 0. The spectral form used The evolutionary development of the magnetic field after 
indicates that the magnetogenesis scenario occured after recombination with no source term is well known. In the 
inflation. absence of velocity fields its adiabatic decrease occurs dur- 

ing the expansion of the universe (e.g. Rees & Reinhardt 
3 2. The solution 1972). The situation looks different if there is shear mo- 

tion. In the considered model the gravitational structure 

For the above postulated categories of the spectrum, the 

begins to form at the recombination time. According to 

expression l|17|) becomes 

the solution l|24|l . the magnitude of the magnetic strength 
g{k } r) — r g(k) r~ 3 , (23) may substantially increase. The question arises whether 
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generation of the large scale field is possible and what is 
the asymptotic form of its spectrum. 
The Fourier mode of the total field is 



Ho(k)l 2 
7 • 10 10 



Hk = [H{Q))k + {5H) k = 
ry s V- 4 /(fc) + ?7 S 2 W)(r- 2 - 4r- 4 ). 



(25) 



Then in calculating of the magnetic energy associated with 
the mode Hk, we restrict ourselves to the leading term in 
r, obtaining 



l^| 2 -^ S 4 f 2 r- 4 |/„(fc)| 2 . 



(26) 



The resulting spectrum evolves through interaction be- 
tween the velocity and magnetic modes which changes its 
initial form. The tendency of these changes is depicted 
by verification of two particular magnetic spectral forms. 
The analytical calculation has been performed for |/„(fc)| 2 , 
(n = 0, 1). The resultant energy density spectra are pre- 
sented below. For a given value of parameters k max , E and 
for kmin = we show in Fig. 1 and 2 the k-dependence of 
\I n (k)\ 2 - the quantity proportional to the magnetic energy 
density of the mode Hk- Two distinct features emerge: an 
energetic discrimination of \I n (k)\ 2 as a function of k and 
the large amplitudes of \H(k)\ 2 oc (3 2 \I n {k)\ 2 ~ (3 2 10 10 . 
Depending on the initial magnetic spectrum, the magnetic 
energy is moved either to smaller and/or larger k. In the 
first case, there is an amplification of the 'magnetic noise' 
by shearing and the energy is carried to small and large 
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Fig. 1. The quantity |io(fc)| 2 proportional to the magnetic 
energy density spectrum as a function of the wave number. 
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Fig. 2. The same as above for |/i(£;)| 2 . Axis of ordinates 
is in log scale. 

which qualitatively shows the results of stretching. In both 
cases the hydrodynamic spectrum boosts the energy to 
different wave intervals. This in turn implies a correlation 
between the magnetic field and the density /velocity fields 
in structures. 

5. Conclusions 



scales. For the initial power spectrum oc f 2 oc k 2 the en- 
ergy transfer goes backward i.e. from larger wave numbers In the presented model magnetic changes result from the 



to smaller wave numbers forming a sort of inverse cascade coupling of velocity shear with the primordial magnetic 
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field. The former appears due to the rapid initialization of Fasc XXV, Les Equations de la Gravitation Einsteinienne, 
structure collapse. The direct implications are (Gauthier-Villars, Paris) 

Grasso, D., & Rubinstein, H., 2001, Phys. Rep., 348, 163 

— field amplification and magnetic energy shift towards 

Hawking, S., & Ellis, G., 1973, The Large Scale Structure 

the larger scales (inverse cascading) and lead to much 

of Space-Time (Cambridge University Press, London) 

steeper magnetic spectra than the initial injection 

Hogan, C, 1983, PRL, 51, 1488 

spectra. 

Kim, E., Olinto, A., & Rosncr, R., 1996, ApJ, 468, 28 

— The shape of the final magnetic spectrum is a func- 

Rees, M., & Reinhardt, M., 1972, A&A, 19, 189 

tion of the initial magnetic spectrum (and implicitly 

Siemicniec, G., & Woszczyna, A., 2004a, A&A, 419, 801 

depends on the primordial density fluctuations). 

Siemieniec, G., & Woszczyna, A., 2004b, A&A, 414, 1 

— The primordial field represents here the indispensable 

Wasserman, I., 1978, ApJ, 224, 337 

seed field for the future amplification (dynamo). Its 
strength at different length scales may differ by several 
orders of magnitude. 

— Also a debate about the top-down versus bottom-up 
scenario for the magnetic field origin in structures may 
be rather unfounded. The seed field may appear simul- 
taneously at all structure scales with strongly diverse 
strength. 

Acknowledgments 

I am grateful to Zdzislaw Golda and to Andrzej Woszczyna 
for discussions and for reading of the manuscript. This 
work is supported in part by 'Komitet Badan Naukowych' 
project 'Cosmic magnetic fields'. 

References 

Berezhiani, Z., & Dolgov, A., 2004, Astropart. Phys., 21, 
59 

Darmois, G., 1927, Memorial de Sciences Mathematiques, 



